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We investigate the transport behavior of polyene molecules sandwiched between two metallic 
contacts using the non-equilibrium Green's function formalism. We calculate both current and 
noise power as a function of applied voltage and show that they decrease with increasing size of the 
polyene molecules. We find that even with symmetric connection to metallic contacts, current verus 
voltage curves can be asymmetric for asymmetrically substituted polyenes. Most importantly, we 
demonstrate a cross-over from Poissonian to sub-Poissonian behavior in the shot noise as a function 
of applied voltage. The algorithm for noise power calculation can be used for designing molecules 
with low noise. 



Recent progress in molecular electronics has led to 
many interesting experiments on single-molecule junc- 
tions [0. In particular, current-voltage (I-V) curves for a 
number of molecules have been measured in lithograph- 
ically fabricated break junctions |^ and using scanning 
tunneling microscopy (STM) |]. Several different the- 
oretical approaches have been suggested for modeling 
charge-transport across single molecules Q, in conjuga- 
tion with molecular modeling at various different levels 
of sophistication, ranging from the simple Hiickel model 
1^ to density functional theories . The nonequilibrium 
Green's Function (NEGF) formalism [Q, which allows 
self-consistent calculations of the electrostatic potential 
along the junction and transmission probabilities, offers a 
very promising method of calculating I-V characteristics 
of single molecule junctions. 

Existing calculations of transport in organic molecu- 
lar junctions have mostly been limited to determinations 
of I-V characteristics. It is well known, however, that 
thermal and quantum fluctuations play important roles 
in the nanometer scale. For future device applications, 
it therefore is essential to investigate current fluctuations 
in these systems. In order to design molecular junctions 
with the lowest possible current fluctuations, an efficient 
algorithm to calculate the noise power in molecular junc- 
tions is needed. To the best of our knowledge, there exists 
curently only one such calculation of the shot noise, in 
Si atomic wires, using a field-theoretic approach |8|. In 
the present Letter we investigate polyene molecules sand- 
wiched between two metallic contacts using the NEGF 
technique. We calculate both the current and the noise 
power as a function of the applied voltage. The latter 
is calculated using an an expression due to Blanter and 
Biittiker |^ . Our I-V curves shows the step- like behavior 
that has been observed in resonant tunneling junctions 
. We show that the current and noise power decreases 
with increasing size of the polyene molecules. We also 
show that I-V characteristics can be asymmetric even 
with symmetric connections to the metal, provided that 
the molecule under investigation is asymmetric One 
of the most interesting results of our calculation is Pois- 
sonian to sub-Poissonian crossover in the shot noise as 



a function of the applied voltage, which we hope can be 
tested experimentally. 

The expressions for the current (/) and the noise power 
S for the general two-terminal device are given by |^ 
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S = \ ( dt{AI{t)AI{0) + AI{0)AI{t)) 

= dE{T{E)[fL{l - h) + fnil - fR)] 

+T{E)[l-T{E)]{fL~fRf}, (2) 

where T{E) is the transmission probability, and 
are the Fermi distributions for left and right terminals, 
respectively, and AI{t) — I{t) — (/). Note the additional 
factor of two due to the summation over spin degrees of 
freedom. 

The principal problem in calculating transmission 
probabilities through a metal-molecule-metal system is 
incorporating the infinitely many degrees of freedom of 
the metallic contacts. The NEGF formalism allows trac- 
ing out these degrees of freedom, such that the dimen- 
sion of the Green's function for the metal-molecule-metal 
system is comparable to the degrees of freedom of the 
molecule [Q. According to the NEGF formalism the 
transmission probability is given by 

TiE)^ trsice[r lGTrG^ (3) 

where G ^ [{E + 0+)I -H T.r]'^ is the Green's 

function, Ti^/j^ = ~ S^^^^J and Sl/_r is the self- 

energy of the left/right terminal. 

The molecular systems we consider are linear polyenes, 
described within the simple Hiickel Hamiltonian, 

H = - CiC^^c,,^ - ^ (4„c,+i,^ + h.c.) (4) 

where c| ^ creates an electron with spin a on site z of a lin- 
ear chain and is the site energy. The nearest neighbor 
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electron hopping integrals correspond to the alter- 

nate double and single bonds in the polyene molecules, 
and are taken to be - 2.6 eV and - 2.2 eV, respectively. 
The site energies are equal for the simple polyene (and 
hence may be taken to be zero). We will simulate asym- 
metry effects (see below) by considering unequal site en- 
ergies, which can be achieved experimentally by chemical 
substituition. Self-energies that account for the semi- 
infinite leads with outgoing plane waves are now calcu- 
lated using the procedure described in Reference . We 
set the hopping integral within the lead - 5 eV, the Fermi 
energy to 5.5 eV and temperature 300 K. This value of 
the Fermi energy corresponds to that of gold electrodes. 
While we report the results of calculations done with this 
one set of parameters only, we emphasize that our calcu- 
lations were done with the lead hopping integral ranging 
from -3.5 eV to - 5.5 eV, Fermi energy ranging from 3.5 
eV - 5.5 eV and temperature IK- 350 K. In all cases 
our results were qualitatively the same. 

Assuming that the voltage drops linearly along the 
polyene molecule, which is a good approximation for 
molecules with small cross section [|o), we have calcu- 
lated the current and noise power for polyenes with 4 
and 10 carbon atoms (butadiene and decapentaene) . In 
Fig. 1(a) we plot current and noise power as a function 
of the bias voltage (VJ,) for butadiene. 

The step-like behavior for the current has been ob- 
served in a number of mesoscopic conductors and is 
a result of resonant tunneling phenomenon Since 
electrons can travel from one terminal to another using 
molecular levels only, a sharp increase in current occurs 
whenever Fermfi level of the metal aligns with the molec- 
ular level. 

Our results for the noise power indicates Poissonian 
{S — 2el) to sub-Poissonian {S < 2el) shot noise 
crossover as a function of the applied voltage. We note 
that sub-Poissonian shot noise has also been observed in 
experiments on mesoscopic conductors 

We introduced asymmetry in the butadiene molecule 
by taking ei = £2 = —0.5 eV and £3 = £4 = 0.5 eV. This 
resulted in asymmetric I-V and noise curves (see Figs. 
1(b)). Note that for positive voltages, noise power ex- 
hibits a nonmonotonic behavior. It would be interesting 
to have experimental confirmations of this result. 

In Reference [|l^ , a different mechanism for asymmet- 
ric I-V curves, based on asymmetric metallic contacts, 
was proposed. A recent experiment has demonstrated 
asymmetric I-V curves for an asymmetric molecule 
Thus, our caculations confirms theoretical findings in 
Reference |l^, that asymmetry in I-V curves can result 
from asymmetry of the molecule too. 

Our results for the polyene with 10 carbon atoms (de- 
capentaene) are shown in Fig. 2. Both the current and 
noise power are now smaller. In addition, many more 
steps now appear in the I-V curve. This can be eas- 
ily explained based on a resonant tunneling phenomena. 
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FIG. 1: Current and noise power for butadiene, (a) symmetric 
molecule, (b) asymmetric molecule. 



Since decapentaene has more energy levels that butadi- 
ene, there are more values of the applied voltage when 
Fermi level of the metal aligns with the molecular levels. 

In Reference ]l3| , where similar I-V curves were calcu- 
lated, it was suggested that low temperatures I-V curves 
can be used for predicting the electronic structure of 
molecules, based on the number of steps in I-V curve. 
Note that the plot of the noise power indicates Poisso- 
nian to sub-Poissonian shot noise crossover similar to one 
we have found for butadiene. 

In summary, we have calculated current and noise 
power for symmetric and asymmetric polyene molecules 
sandwiched between two metallic contacts using the 
NEGF formalism. Our results indicate that the cur- 
rent and noise power decrease with increasing size of 
polyene molecules. We have shown that asymmetry in 
I-V and noise power curves can result from asymmetry 
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FIG. 2: Current and noise power for polyene with 10 carbon 
atoms. 



of the molecule. We have also shown Poissonian to sub- 
Possonian shot noise crossover as a function of the ap- 
plied voltage. With constantly shrinking size of electrical 
circuits, the limit of single molecular junctions is being 
rapidly approached. For single-molecule junctions exper- 
minctal results on the noise power are as important as I-V 
characyeristics. We believe that the algorithm presented 
here for the noise power calculation can provide an ef- 
ficient theoretical approach to designing molecules with 
the lowest possible noise. 

This work was supported by NSF DMR-0101659 and 



[1] 
[2] 
[3] 

[4] 

[5] 
[6] 
[7] 

[8] 
[9] 

[10] 

[11] 
[12] 
[13] 



M.C. Petty, M.R. Bryce and D. Bloor (eds.). Introduction 
to Molecular Electronics (Oxford University Press, New 
York, 1995). 

J. Riechert, R. Ochs, D. Beckmann, H.B. Weber, M. 
Mayor and H.v. Lohneysen, Phys. Rev. Lett. 88, 176804 
(2002). 

V.V. Shorokhov, P. Johansson and E.S. Soldatov, J. 
Appl. Phys. 91, 3049 (2002); F. Moresco, G. Meyer, K.- 

H. Rieder, H. Tang, A. Gourdon and C. Joachim, Phys. 
Rev. Lett. 86, 672 (2001); 

J. Guo, M. Lundstrom and S. Datta, Appl. Phys. Lett. 
80, 3192 (2002); C. Krzeminski, C. Delerue, G. Allan, D. 
Vuillaume and R.M. Metzger, Phys. Rev. B 64, 085405 
(2001); M. Paulsson and S. Stafstrom, Phys. Rev. B 64, 
035416 (2001). 

I. R. Peterson, D. Vuillaume and R.M. Metzger, J. Phys. 
Chem. A 105, 4702 (2001). 

M. Di Ventra, N.D. Lang and S.T. Pantelides, Chem. 
Phys. 281, 189 (2002). 

S. Datta, Electronic Transport in Mesoscopic Systems 

(Cambridge University Press, UK, 1995). 

Y.-C. Chen and M. Di Ventra, |cond-mat/0208'066| . 

Ya. M. Planter and M.Biittiker, Phys. Rep. 336, 1 

(2000). 

A. Nitzan. M. G alperin. G.-L. Ingold and H. Grabert, 



physics/0207124 



E.G. Emberly and G. Kirczenow, Phys. Rev. B 64, 
235412 (2001). 

J. Heurich, J.C. Cuevas, W. Wenzel and G. Schon Phys. 

Rev. Lett. 88, 2568 03 (2002). 

N.A. Zimbovskaya, |Eond-mat/0208304 



